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Abstract 

Exact critical properties of the one-dimensional SU(iV) interacting fermion 
model with open boundaries are studied by using the Bethe ansatz method. 
We derive the surface critical exponents of various correlation functions using 
boundary conformal field theory. They are classified into two types, i.e. the 
exponents for the chiral SU(iV) Tomonaga-Luttinger liquid and those related 
to the orthogonality catastrophe. We discuss a possible application of the 
results to the photoemission (absorption) in the edge state of the fractional 
quantum Hall effect. 
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I. INTRODUCTION 



One-dimensional (ID) quantum many-body systems with open boundaries have attracted 
much interest recently in connection with various problems in condensed matter physics such 
as the Kondo problem, the tunneling in a quantum wire and the edge state of the fractional 
quantum Hall effect (FQHE). A number of exactly solvable models with open boundaries 
have been known so far, e.g. the XXZ Heisenberg model JIJ, the interacting boson model 
0, the interacting fermion model ||, the Hubbard model |fj], and the 1/r 2 quantum model 
||, etc. Quantum impurity models such as the Kondo model and the Anderson model also 
have a deep connection with the boundary problem |||7]]. All the above systems exhibit the 
surface critical behavior near the boundary. The corresponding surface exponents, which 
should be different from the bulk ones, can be obtained by using the finite-size scaling in 
boundary conformal field theory ||. 

In this paper, we obtain the exact surface critical exponents for the ID SU(iV) interacting 
fermion model with open boundaries by combining the Bethe ansatz solution and boundary 
conformal field theory. By examining boundary conditions carefully, we classify the surface 
exponents into two categories, i.e. the exponents for the chiral Tomonaga-Luttinger liquid 
and those related to the orthogonality catastrophe. The latter is shown to be related to the 
X-ray problem in ID chiral systems. We apply our results to the edge states of the FQHE, 
and predict some expected behaviors for the photoemission (absorption) singularity. Some 
of the results obtained in this paper were previously conjectured in J7|. 

The plan of the paper is as follows. In sec. II, we introduce the continuum fermion 
model with SU(iV) spin symmetry in open boundary conditions, and give the Bethe ansatz 
equations. In sec. Ill, based on boundary conformal field theory, we derive critical exponents 
for various correlation functions from the finite-size spectrum computed by the Bethe ansatz 
solution. In sec. IV, we discuss a possible application of the present results to the X-ray 
photoemission (absorption) problem in ID chiral electron systems. A brief summary is 
given in sec.V. 
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II. MODEL AND BETHE ANSATZ SOLUTION 



We consider the interacting fermion model with SU(iV) spin symmetry in open boundary 
conditions. The mutual electron interaction is of 5-function type, and boundary potentials 
are introduced at both ends of the open chain. The Hamiltonian is thus given by 

Q2 N N m 



N N m 



2c y: s( 

i<j,m,n 



XXwH^m) + lL^{Xi,m ~ L)), (C > 0) 



m=l i=l 
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where N m is the number of electrons with spin index m (= 1, 2, • • • , N) of SU(iV) symmetry 
and L is the system size. The last two terms represent boundary potentials with coupling 
constants 7o(l). We note that this type of the SU(iV) fermion model with 5-function interac- 
tion was first solved by Sutherland under the periodic boundary condition many years ago 
||. We now wish to solve the above model ([J) for the open chain with boundary potentials 
7o(l)- This can be preformed by generalizing Gaudin's method developed for the boson 
model with open boundary conditions 0. The main idea is that one can treat the open 
boundary problem more easily by introducing fictitious mirror-image particles with respect 
to the boundary. We then diagonalize the Hamiltonian (JJ) following standard Bethe ansatz 
techniques developed for periodic systems M, and end up with the basic algebraic equations 
for rapidities kj and A£>, 
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(p = 0, L) are the phase shifts due to the boundary scattering. Quantum numbers Ij and 



J$ are positive integers (or half-odd integers) which classify the elementary excitations. 
The total energy is simply given by 



N 
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Note that the above equations contain the terms with arguments like A^ + A^, reflecting the 
introduction of image particles @-|J], in contrast to the ordinary Bethe ansatz equations for 
the periodic case. These equations for the SU(2) case were obtained firstly by Woynarovich 
||. Defining new variables, A_ Q = — A^, we can rewrite the above equations into more 
tractable form, 

-2(%-A«)- 



/2k \ Ml 

2Lkj + <po(kj) + (p L (kj) = 2-Klj + 2 tan" 1 ( — ) - J2 2 tan_1 

^ ' P=-Mi 



2 tan" 1 

Mi 

+ 2tan 

P=-M t 

Mjv_ 2 

53 2 tan" 1 

{3=-M N _2 



2(A«-Ag" 1) ) 



2tx J n — 2 tan 



2 tan" 



c 

i/2A(0x 



(5) 



rA« - A<? 



2(A(f - 1 ) - AjJ 



c 

(N-2), 



l+l 



2 tan" 

P=-M l+1 

27r J n — 2 tan 



v c y 1 \ c j 

2(AW - Ag +1 ^ 



1 < / < AT - 2, (6) 



+ 53 2 tan" 1 



/3 



(7) 



After this transformation, the structure of the equations formally resembles to that for the 
periodic case ||. This fact makes the following analysis much easier . 

In the following, we will be concerned with the case of repulsive boundary potentials 
(7 P < 0) for simplicity. All the rapidities then turn out to be real in this case. Taking 
the thermodynamic limit, we now recast the algebraic equations fl5|) ~(0) into the integral 
equations for the density functions of rapidities, 
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where A+^ = — A_ = kp\ which is determined by the condition, 
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In the absence of magnetic fields, one can see that Ap = -A W -> +oo (1 < / < N- 1), and 
the system recovers SU(A) spin symmetry. It is convenient to divide the density functions 
into bulk and boundary parts, 

p(k) = Pbuik(k) + —pb(k). (12) 

In zero magnetic field, the coupled integral equations for the density functions can be reduced 
to simple ones for p bu i k and p b from eqs.(|) ~(|10D by using Fourier transformation. The 
results are 
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The total energy is given by 
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Eqs.([lq) ~ flTTD determine the energy spectrum of the model. It should be noticed that in 
the thermodynamic limit, the bulk properties of the present model is the same as those for 



the periodic model ||, which have been already studied in detail [ I0| . To examine critical 
properties, we thus need to study the finite-size spectrum, which will be done below. 



III. BOUNDARY CRITICAL PROPERTIES 

In this section, we compute the energy spectrum for the finite system, and discuss low- 
energy critical properties using boundary conformal field theory. We then obtain the surface 
critical exponents of various correlation functions. 

A. finite-size spectrum and conformal properties 



Applying a standard technique in the Bethe ansatz method pj| , we now obtain the 
finite-size energy spectrum from the basic equations derived in the previous section. For 
this purpose, let us first express the total energy, eq.([Tj]), in terms of the dressed energies, 
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where the dressed energies e"> are determined by the following integral equations, 
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The equivalence of the two expressions (|T7|) and (p~8f ) can be easily checked by directly 
comparing them after formally solving the integral equations by the iteration scheme. The 
formula (|18|) is particularly useful to compute the excitation spectrum. 

Let us start with the corrections to the ground state energy. By directly applying the 
Euler-Maclaurin formula, 

1 & fn\ r(P+i/2)/N i / ,/ a -l/2\ ,/6 + l/2\\ 
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to eq.([|), we easily find the finite-size corrections to the ground state energy p2| , 
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which is scaled by the velocities t>/ (Z = 0, 1, • • • , N — 1) defined by 
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with a(°\A^) = By exploiting the finite-size scaling for open boundaries fL3|,Sf, we 

can see from eq.(p3|) that the Virasoro central charge for the charge sector (/ = 0) is given 
by c = 1. On the other hand, all the velocities of the spin excitation takes the same value 
(Z = 1, 2, • • • , N — 1) in the absence of magnetic fields, and thus the central charge for the 
spin sector turns out to be c = N — 1, namely, the rank of SU(iV) Lie algebra. We shall 
study conformal properties in more detail by examining the excitation spectrum below. 

Elementary excitations are classified into two types, i.e. those for primary fields and 
for descendant fields. The former can be described by the excitations which change the 
number of particles. They are computed by changing the cut-off parameters in ([Tj|) as A± — > 



A± + AA±*. On the other hand, the excitations for descendant fields are simply given by the 
particle-hole type excitations with the fixed number of particles. These manipulations are 
performed straightforwardly, and we end up with the finite-size spectrum for the excitation 
energy, 
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with V = diag(t>o, v i, fjv-i), where n+ are non-negative integers denoting particle-hole 



excitations. Here £ is the N x N matrix of the so-called dressed charge ||14j| , whose compo- 
nents £ij = ^(A^) (0 < i, j < N — 1) are determined by the following integral equations, 
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In eq.(^), the quantum numbers classifying elementary excitations are defined as 

AM (!) = AMjp - ^(N-l), 1<1<N-1, 
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Here n& is the number of particles localized at boundaries, which is given by 



i r*P 



(30) 



and AAT^ is an integer specifying charge excitations, whereas AM^'s are integers which 
label (N — 1) kinds of spin excitations. 



S 



We can now read off conformal weights A& from eq.(^), using finite-size scaling ar- 
guments |15|. Surface critical properties near the boundary are determined by boundary 
scaling operators <fi(t). The critical exponent for correlation functions of a boundary opera- 
tor, ((f) (t) 0(0)) ~ l/t x , is given by 

AT-l 
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where the N x N matrix Cf = (£ ) £ is given in the absence of magnetic fields, 



1 NKp ~ N 

-1 
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Here K p = C,o {kp S) )/N is the dimensionless coupling constant for the charge sector (so- 
called Tomonaga-Luttinger parameter), where the dressed charge £oo(^) is determined by 
the integral equation, 



£oo(fc) = l+ / ' dkG(k-k')C 00 {k'). 



(33) 



Note that the system is now regarded to be chiral due to the presence of open boundaries, 
and quantum numbers carrying currents do not appear in the conformal weights ([31]). This 
implies that an effective theory of the present model is given by the holomorphic piece of 
conformal field theory. From eqs.([H|) and (|32"D, we can see that critical properties of the 
charge and spin sectors are respectively described by the U(l) Gaussian model with the 
central charge c = 1 and the level-1 SU(iV) Wess-Zumino-Witten model with c = N — 1. 

The expression ([H]) with (|32"D is one of the main results in this paper. We wish to 
emphasize that this formula is quite general and is applicable to many other SU(iV) quantum 
models with boundaries such as the t- J model, the Hubbard model 01, etc. 
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B. surface critical exponents 



We are now ready to obtain the surface critical exponents of various correlation functions. 
As is seen from eq.(^), the effects of boundary potentials are just to shift the quantum 
number as AM"' — > AM"' — jt{N — I). One readily notices that such an effect of boundary 
potentials is equivalent to that of twisting boundary conditions, which does not change the 
critical exponents in general. Thus when we determine the critical exponents from (|5ID, 
we should discard n^-dependence in AM by redefining the quantum number as AM^ — 
j&(N — I) — > AM®. Therefore, for the long-time behavior of the single-particle Green 
function, (c£(t)c o (0)) ~ l/t n , we obtain its critical exponent by setting ANh = 1, AM® = 
(1 < I < N- 1), 

„ = ^ + ^I. (34) 
' NK p N K J 

Furthermore it is seen that the density- density correlation function and the spin-spin cor- 
relation function show the following asymptotic behavior, by taking ANh — 0, AM|' = 



[1 < I < N- 1), and n 



(O(t)O(O)) ~ constant + 4- (35) 

t 2 



Namely, this asymptotic behavior is determined by descendants of the primary field with 
A5 = 0, and hence the critical exponent takes the canonical (integer) values. The anomalous 
exponent appears only in the single-particle Green function. These characteristic proper- 
ties are inherent in the chiral Tomonaga-Luttinger liquid |H ] , which is quite contrasted to 



ordinary periodic systems as will be seen in the next section. 

We wish to stress that the formula (|31~1) for conformal dimensions possesses another 
important information for boundary critical properties related to the orthogonality catas- 
trophe. This is realized when one considers the problem in which the boundary potentials 
are time-dependent |22|,|7|. For example, suppose that the boundary potentials are sud- 
denly turned on at to- Then the long-time asymptotic behavior of correlation functions 
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(0(t2)0(ti)) for ti <C to *C £2 should show quite different properties from (|3~4f) and ([35]). In 
this case, the phase shift caused by boundary potentials plays an essential role, and then 
n& cannot be ignored by redefining the quantum numbers. Therefore, retaining n& in (|5ID , 
we get anomalous exponents for various correlation functions. We show several examples 
below. 

(i) single-particle Green function: (4(t)c Q (0)) ~ 1/f. Putting AN h = 1, AM? = 
(1 < I < N- 1), we have 

1 , ^2 N-l , , 

" = NK P { - 1 - *>' + — (36) 

(ii) density- density correlation function: (n(t)ra(O)) ~ l/t Qc , = Z)a c l c a)- Taking ANh 
I). AMf = (1 < / < JV - 1), we then obtain 

°-=m;- (37) 

(iii) spin-spin correlation function: (S a (t)S a (0)) ~ (£" = J2 a ,/3 c a T a/3 c /3 with being 
a fundamental representation of SU(iV)). Putting AN h = 0, AMf = (1 < / < iV - 1), 
we have 

a, = (38) 

Thus the surface critical exponents depend not only on the Tomonaga-Luttinger param- 
eter K p but also on the fractional number n& of localized particles at the boundary. We will 
see in the next section that these critical properties have a deep connection to the X-ray 
edge problem in ID chiral systems. 

We can evaluate the above exponents easily in some limiting cases. In the case of 
noninteracting electrons (c = 0), we see K p = 1 from ([33]), and the critical exponents (|34|). 
(|35|) (|36D, fl3~T|) and (|38| ) are all reduced to those obtained previously for the single- impurity 
SU(iV) Anderson model with infinitely strong Coulomb interaction 0. On the other hand, 
in the strong-coupling limit c — > +oo or low-density limit kp — > 0, it is easily seen from 



eq. (p3[) that K p = 1/N, and so the exponent of the single-particle Green function eq.(p4| 
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is 77 = 2 — 1/N. In general, K p may change in the range [1,2 — 1/N] according to the 
interaction strength and the density of particles. 



C. comparison with the periodic model 

To conclude this section, we compare the above results with bulk critical exponents under 
periodic boundary conditions. The finite-size spectrum for the SU(iV) model in the periodic 
boundary conditions is given by [|L4] 

I = x(l AMT ^ 1 ) Ty ^ 1 ) AM + DT ^ TD + ^ 1 ^^ ) +^0' ( 39 ) 

where the newly introduced quantum numbers, D T = (D , D 1; .D^-i), are given by 

D = 1 , (mod 1), (40) 



A = \ , (modi), Z = l,...,iV-l, (41) 



with AM^' = 0. These quantum numbers carry the large momentum transfer JI~4 



Using finite-size scaling arguments for periodic systems p5| , we have the conformal 

dimensions of scaling operators, 

1 N-l 
A + + A_ = -AM T C / AM + D^D + ^ (n? + n (l) ). (42) 



4 

From this formula, we obtain critical exponents. We list several examples in the following, 
(i) The single-particle Green function: 

1/1 „\ iV-1 



p 



' = SfE + *' + — ■ (43) 



In the strong-coupling limit or low-density limit, this reduces to [ PH , [H| ] 

3 1 1 

^=2-iV + 2A^- (44) 
(ii) The 2A r /ci?-oscillating term (kp: Fermi momentum) of the density-density correlation 
function: 
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a c = 2NK P . (45) 



We have a c = 2 in the strong-coupling limit or low-density limit, 
(hi) The 2/cp-oscillating term of the spin-spin correlation function: 

K p + N-1 



which reduces to 



<*. = , (46) 



a - = 2 ^-h + w^ (47) 



in the strong-coupling limit or low-density limit [|16[ 

Comparing these results with the previous ones, we can see that the bulk critical expo- 
nents are quite different from the boundary ones, though bulk thermodynamic properties 
should exhibit the same behavior in both cases. 

IV. APPLICATION TO THE EDGE STATE OF FQHE 

In this section, we briefly discuss a possible application of the results obtained in the 
previous section to the edge state of the FQHE as a typical example of ID chiral fermion 



systems [|nj] . We first note that if an appropriate value is chosen for the Tomonaga-Luttinger 
parameter K p , the above SU(iV) model with open boundaries can give an effective theory for 
the edge state of a certain hierarchy in the FQHE. Namely, the edge state of the fractional 
quantum Hall effect with filling v = N/(Nm + 1) (m even) can be modeled by the above 
open system by choosing K p as K p = u/N J7|. In fact, one can easily check that by this 



choice of K p , the formulas (p3) and (B5|) exactly reproduce the critical exponents for the 



chiral liquids proposed for the edge states [^1[ . 

We shall focus on the X-ray photoemission (absorption) problem in edge states. The 



Fermi-edge singularity problem in ID electron systems has attracted current interest flq-pO 
In the edge state of the FQHE, electrons move only in one direction and the backward scat- 
tering due to impurities is irrelevant, and hence the system can be treated as a chiral system. 
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In experiments of the X-ray photoemission or absorption, the core hole may be suddenly 



created, resulting in a problem with time-dependent boundary conditions [|7|,|22| . That is, for 
t < to the boundary is free, and at t = t the boundary potential suddenly switches on. Then 
bulk electrons show critical low-energy behavior inherent in the orthogonality catastrophe. 

In the X-ray absorption problem, one electron is excited in the final state. Thus putting 
ANh = 1, AM^ = (1 < / < TV — 1), we have the critical exponent for the X-ray absorption 
in this system, 



1/ NS\ 2 N-l 

a *=-A 1 --) + — < 48) 

for the filling v = N/ (Nm + 1) with even m, where S is the phase shift caused by the localized 
electrons which screen the core hole potential. On the other hand, the critical exponent for 
the photoemission is obtained by taking ANh = and AMf = in eq.([H|), because one 
hole carrying neither charge nor spin is generated in the final state, 

°* = -{-)- (49) 

We expect that such anomalous exponents may be observed in the X-ray photoemission or 
absorption experiments for the edge state of the FQHE of the filling v = N/ (Nm + 1) with 
even m. 



V. SUMMARY 

We have studied exact boundary critical properties of the SU(iV) interacting fermion 
model with open boundaries by using the Bethe ansatz solution and boundary conformal 
field theory. It has been shown that the surface exponents which govern the critical behavior 
near the boundary depend on the two continuously varying quantities, i.e. the dimensionless 
Tomonaga-Luttinger parameter K p and the fractional number of localized electrons at the 
boundary. We have also discussed a possible application of the results to the Fermi-edge 
singularity problem in the edge state of the FQHE as a typical example of ID chiral systems. 
The exact exponents for the X-ray absorption and photoemission have been derived. 
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